The shape of the nucleon 
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We show that all four of the nucleon form factors can be very well explained using the manifestly 
covariant spectator theory. The nucleon is modeled as a spherical state of three constituent quarks 
with electromagnetic form factors, all determined by the fit to the data. 



There have been many theoretical attempts to explain 
the recent high precision polarization transfer measure- 
ments of the ratio of proton form factors Ge^/Gmp [HQ 
and neutron form factors Geu/Gmu [3 Q As shown 
in Ref. ||, recent calculations that model the micro- 
scopic structure of the nucleon using the Skyrme model 
, front-form quantum mechanics 0, 0, E3 1 or point- 
form quantum mechanics [Tl| fail to fit all of the form 
factor data. Models inspired by vector dominance with 
7 to 12 parameters do much better but do not 

model the quark structure of the nucleon. These difficul- 
ties have lead some to suggest 0, Q that the proton is 
deformed, and this possibility has been widely discussed, 
reaching Science News [HI and USA Today While it 
is possible for a spin 1/2 system to be deformed 0], the 
accompanying rotational bands and large spectroscopic 
factors have not been observed. 

In this letter we use a simple model to show that the 
beautiful Jefferson Lab data do not require that the pro- 
ton be anything other than spherical. Our model assumes 
that the nucleon is composed of three valence constituent 
quarks (CQ), and that the main role of the surround- 
ing sea of gluons and quark-antiquark pairs is to dress 
these CQ. Through this dressing the CQ acquire mass, 
size, and structure, and are composite systems with elec- 
tromagnetic form factors of their own. The structure 
of these form factors is given by QCD, which requires 
that the quarks become point particles as the momentum 
transfer Q 2 — > oo, and that the 7* — > qq spectral function 
at finite Q 2 be the sum of the qq bound state poles and 
resonances allowed by confinement (vector dominance). 

We use the covariant spectator formalism 

El 

to model 

the bound states of three CQ. Using this method, all 
of the Poincare transformations (including rotations and 
boosts) can be handled exactly, important when dis- 
cussing the shapes of states. This formalism has already 
been used successfully to calculate the relativistic wave 
function of 3 H |l8|, UjJ |2fJ , which is described by a vertex 
function with two of the three hadrons on mass-shell, and 
the third off-shell. The wave function has the form 



with m q the mass of the (dressed) quark, P = p\+p 2 +p 3l 



(j>2,P3\T\P) 
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r the vertex function, and p 2 



pi 



We adopt a 



simple approximation in which the distributed mass of 
the spectator diquark system, p 2 = (p2 + P3) 2 , is fixed 
at a mean value of m 2 . Both m s and m q disappear from 
the final result. 

If the mass of the nucleon is larger than the mass of 
the three quarks, M > 3m q , the pole in the propagator of 
the off-shell particle 1 moves to the real pio axis. If this 
pole were not cancelled by a zero in the vertex function, 
it would describe free scattering of the quarks. How- 
ever, it has been shown that a dynamical consequence 
of confinement is to automatically produce a zero in the 
vertex function that cancels the pole [2l| . In this work 
we include confinement by modeling the wave function di- 
rectly, assuming that Q is smooth at the point p\ = m 2 . 

In the symmetric quark model the wave function can 
be written as the sum of two components 



NR 



+ sin Oififi) V, 



(2) 



where ip is a scalar function symmetric in the coordinates, 
and / = ±5, s = ±| are isospin and spin projections, 
respectively. The states 4>j couple quarks 2 and 3 to total 
isospin J = or 1: 
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with xl = "l, Xi 2 = di, £° 3 = -j= (u 2 d 3 - d 2 u 3 ), and 
C 3 + = u 2-" 3 , £33° = ^ (u 2 d 3 + d 2 u 3 ), and ^~ = d 2 d 3 . 
The spin wave functions have the same form. The sum 
P|) is symmetric under the interchange of quarks 2 and 
3, and when cos 9 = sind, Eq. J2J) is also symmetric un- 
der interchange of all the quarks, so in the following we 
discard the quark labels 1,2, and 3. 

Following J5J), we assume that the relativistic wave 
function Q has the manifestly covariant form 

+ ^=(757p)»„^^(^s)(r.e 1 )x 7 V'i( J P,p) (4) 
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where p is the on-shell diquark four-momentum (p = 
nig), u(P, s) is the spinor of the nucleon, a is the Dirac 
index of the off-shell quark, and rj is the covariant polar- 
ization vector of the spin-one diquark. The scalar func- 
tions ijio and tp\ can depend on only one variable, which 
we choose to be 



M 2 



X 



(P-P) 2 



Mm, 



(5) 



In the cm frame, \ — 2E s (k), where E s (k) = yl + n 2 
and k are the energy and magnitude of the three- 
momentum p of the on-shell diquark in units of m s . Note 
that the wave function depends only on the magnitude of 
p, and hence is spherically symmetric. To get the correct 
asymptotic dependence of the nucleon form factors 
predicted by QCD, we choose a simple double pole de- 
pendence for the scalar function ipo 



MPp) = 



J (/3 1 -2 + X )(^2-2 + x) 



(6) 



where 0i and 02 are range parameters in units of Mm B 
and No is a normalization constant. In the nonrelativistic 
(NR) limit (k « 1), ipo becomes the familiar Hulthcn 
wave function 



N 



nr m s (pi + K 2 )(f3 2 + K 2 ) 



(7) 



We chose the form of ipi to be related to ipo by a rela- 
tivistic factor that approaches unity in the nonrelativistic 
limit 
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This choice insures that both the proton and the neu- 
tron charges are correctly given by a single normalization 
condition (see Eq. (|21(l below), and is a relativistic gen- 
eralization of the symmetry requirement ip\ = ipo- With 
these definitions, the NR limit of ^ a in its cm frame 
becomes $nr- 

With these ansatze, the form factors of a nucleon with 
isospin / are given in the spectator theory by the rela- 
tivistic impulse approximation (RIA) [22T | 



J? = Fi(Q 2 ) + F2(Q 2 ) 



2M 

= 3u{P+,\')Jfu(P-,X), 



(9) 



where F{ are the nucleon Dirac and Pauli form factors 
and r j is the current of a single quark with isospin / 
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7,75^^757,' W V MP+,P)MP-,P) , (10) 



and we define P± = P ± ^q. The CQ current, jj , is 
defined in Eq. below. We have summed over the 
polarization vectors of the diquark, using 



■q rj — A = 



(11) 



The total result is three times the contribution of one 
quark. 

The manifestly covariant current is most easily 
calculated in the Breit frame with Qo = Q/M, Pq = 
y/1 + Ql/4 and P± = M{P , 0, 0, ±Q /2}. In this frame 
■00 becomes 

N 



^Po t = MP±,p) = 



(12) 



m s D^{Q) 

with D±{Q) = (/?! - 2 + x ±)(Pa - 2 + X ±) and X ± = 
2PoE s (k) ± kQqz, where z is the cosine of the angle 
between p and q. Similarly, the ipt wave functions are 
obtained from Eq. (|SJ) by substituting \ ~ * Xt ■ 
The most general form of the CQ current is 



Ji 
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where j\ and j% both depend on isospin (T3 = ±) and Q 2 
through the CQ form factors 



ii = /i+ + T3/1- ; 
h = f%+ + 7-3/2- ; 



/l± — T \fld 

fa± = 5 (M«/2« ± Vdlid) ■ (14) 



Here \i u and /id are the anomalous moments of the u and 
d quarks, adjusted to fit the nucleon anomalous moments. 

Three principals guide the choice of CQ form factors. 
First, QCD requires that the CQ should become elemen- 
tary point-like quarks as Q 2 — > 00, and this requires 
h±{Q 2 ) -> and h±{Q 2 ) -> A± as Q 2 -> 00. To 
assure that the charge assignments remain in the ratio 
of 2/3 to —1/3 at both large Q 2 (when the quarks are 
point-like and the charge assignment is given by QCD) 
and small Q 2 (where the charge assignments are given by 
quark model phenomenology) we require A+ = A_ = A. 
[The charges could then be changed to 2/3 and —1/3 
at Q 2 — > 00 by absorbing A into Nq.] Second, guided 
by vector dominance, we expect the individual isoscalar 
and isovector form factors to have a monopole 
contribution 

G(A2) - J^rm < 15 > 

with 4 scale parameters, Af ± , independently adjusted 
to fit the data. We find the effective vector dominance 
masses to be strongly channel dependent (see Table [fl . 
Third, we allow for direct coupling of the photon to pi- 
ons inside the nucleon by adding a purely isovector dipole 
charge contribution to /i_ 
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Assembling these pieces gives the following 

/ 1+ (Q 2 ) = e+ [(l-A)G(A 1+ ) + A] 

h-{Q 2 ) = e_[(l - A - A 7r )G(A 1 _) + AttGtt + A] , 

/ 2± (Q 2 ) = M±G(A 2± ) (17) 

where e + = i, e _ = i, /i± = ± Md)- 

To complete the calculation the second term in <|10f) 
must be reduced. The isospin algebra yields 



J{i+2) 



(18) 



where i — {1,2}. Averaging over the asimuthal angle 
and using the Dirac equation allows the extraction of the 
nucleon form factors F\ and F 2 
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where R\, R 2 , and -R 3 are 
fl 
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At Q 2 — the relativistic factor © cancels giving 



with D±(0) = D{0) and 
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The condition I\ = 1 now insures that the correct pro- 
ton and neutron charges emerge automatically. In a full 
theory this feature would emerge without additional as- 
sumptions |22t |23| so we do not regard Nq as an ad- 
justable parameter. 

As Q 2 — > oo, the form factors f\± approach Xe± and 
the form factors f 2 ± approach zero as A 2 ±/Q 2 . Then 
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where m± = /i±A 2 ^ and J 2 is 
d 3 K 1 



(2tt) 3 2E s (k) [1 + k 2 (1-z 2 )] 2 



- r 3 e_] I 2 

r 3 m_] J 2 (23) 
0.0179. (24) 
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FIG. 1: (Color on line) The quark form factors. Solid lines are u 
and dashed lines are d form factors; upper 4 lines are /i (Model I 
larger than II) and lower 4 lines /2 (Model I again larger than II). 
Asymptotic values of the /i form factors are shown. 



TABLE I: Ten adjustable parameters (all dimensionless) . The 
11th adjustable parameter, A, is 0.230 (I) or 0.177 (II). 



Model f3i, §2 fiu, Af+, At 



A^ 

A 2+i A 2- 



Ati-, A,. 



(20) 


I 


0.057 


1.125 


7.69 


0.362 


0.245 




1.982 


-0.837 


10.67 


1.82 


2.15 




II 


0.069 


1.126 


3.23 


0.386 


0.244 


iving 


















1.543 


-0.825 


5.81 


1.23 


0.742 



We present two fits to the nucleon form factors. The 
parameters of each fit are shown in Table ^ with prop- 
erties displayed in Table [n] quark form factors shown in 
Fig. ^ an d nucleon form factors shown in Fig. |2 (We 
have not used data for Ge p obtained from Roscnbluth 
separations because it is now believed that such data 
are contaminated by two photon exchange contributions 
\2(fy.) Model I was obtained by minimizing % 2 , and gives 
excellent agreement (see Table ITT)l . Model II starts from 
Model I and further adjusts the parameters to give an 
excellent fit to the charge radii. 

The two magnetic form factors Gm p and Gmu have a 
very similar shape, which seems to disagree with the high 
Q 2 trend for Gmu- However, this "trend" is defined by 
only two data points with large errors. These models pre- 
dict the shape that new, more accurate measurements of 
Gmu (being analyzed |2jJ) should follow. The fits to the 
electric form factors are both excellent, and predict that 



TABLE II: Properties of the two fits. For comparison, the 
experimental values are r 2 , = 0.780(25) and r 2 = —0.113(7). 



Model x 2 /datum r 2 (fm 2 



r 2 (fm 2 



vo 



I 

II 



1.18 
1.94 



0.800 
0.771 



-0.065 
-0.106 



36.39 
60.87 




s 0.1 





FIG. 2: (Color on line) Data for the nucleon form factors compared 
with Model I (solid lines) and Model II (dashed lines). The Gmp 
and G\ fn data were used by Bosted in his global fits to the form 
factors l24l . The Ge p data are from JLab Hall A by Jones, et ai.0] 
(squar es) and Gayou, et (triangles). The Geu data are from 

Bates l25l (circle), and from JLab Hall C by Zhu, et.al. (trian- 
gle), Warren, et (squares), and Madey, et al. ,5J (diamonds). 

Gep will go through zero at Q 2 - 8 GeV 2 . The quality 
of our fits is comparable to the recent vector dominance 
model of Ref. [ij. 

The observation that QF2 P /Fi p is nearly flat for the 
highest points measured has lead to speculation that high 
energy ep elastic scattering violates the conservation of 
helicity |28j, and that the proton may have an exotic 
electromagnetic structure. Fig. [3] shows that our Models 
predict that the flatness currently observed is only tem- 




g(GeV) 

FIG. 3: (Color on line) Double log plots showing the data for 
Q F2p / (p-p F\ v ) and Q 2 F^ v / (fi p Fi p ) versus Q and Q 2 respectively 
(labeled as in Fig.0. The dotted lines are the asymptotic predic- 
tion for Model I, Q 2 F 2p / (u p F lp ) = 7.06, obtained from Eq. l23t . 



porary, and that QF2 P /F\ p does indeed approach zero as 
1/Q, as expected. The asymptotic behavior is delayed 
until Q 2 > 100 GeV 2 . This conclusion is largely inde- 
pendent of the running of the QCD coupling constant 
a s , which we have not discussed. 

We conclude that the new form factor data do not 
provide evidence that the nucleon is anything other than 
a spherically symmetric state of three relativistic con- 
stituent quarks |29|. We also found that the data can 
not be fit unless it is assumed that the CQ have an elec- 
tromagnetic structure. The precision of the data allow 
this structure to be determined empirically (subject to 
uncertainties as shown in Fig.^l. This fit requires a pion 
"cloud" (the term) and leads to the conclusion that 
the onset of perturbative QCD is delayed until large Q 2 , 
as expected. 
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